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We derive a Kondo Lattice model with a correlated conduction band from a two-band Hubbard Hamiltonian. 
This mapping allows us to describe the emergence of a robust pairing mechanism in a model that only contains 
repulsive interactions. The mechanism is due to topological confinement and results from the interplay be- 
tween antiferromagnetism and delocalization. By using Density-Matrix-Renormalization-Group (DMRG) we 
demonstrate that this mechanism leads to dominant superconducting correlations in a ID-system. 
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The origin of unconventional superconductivity remains as 
one the most important open problems of condensed matter 
physics. Physicists do not agree on the mechanism that pairs 
the electrons to form the superconducting condensate. To a 
very good approximation, electrons only interact via the re- 
pulsive Coulomb interaction. Consequently, as it was pointed 
out recently by Anderson the crucial question is: "How 
can this repulsion between electrons be eliminated in favor 
of electron pair binding?" The problem becomes even more 
puzzling if we consider that the Coulomb interaction is big- 
ger than the bandwidth for most of the unconventional super- 
conductors. Then, the first challenge is to demonstrate that a 
"pairing force" can exist in model Hamiltonians that only con- 
tain strongly repulsive interactions. Although different pairing 
mechanisms have been proposed over the last twenty years, it 
is not always clear if those mechanisms actually work or if 
they are robust under the presence of long range Coulomb in- 
teractions. This is mainly due to the lack of controlled approx- 
imations for solving models of strongly interacting electrons 
in 2D or 3D. 

Another aspect that is quite ubiquitous in unconventional 
superconductors is the proximity of the superconducting state 
to an antiferromagnetic (AFM) phase. This observation sug- 
gests that antifeiTomagnetic correlations are related to the 
pairing mechanism. However, although several "magnetic" 
pairing mechanisms have been proposed |2 |, it is still unclear 
how the interplay between AFM correlations and itineracy 
leads to a "glue" that is strong enough to hold the two elec- 
trons together. It is the purpose of this Letter to show how a 
robust pairing mechanism emerges out of this interplay and to 
demonstrate that it leads to dominant superconducting corre- 
lations in a two-band Hubbard chain. Moreover, we will see 
that the pairing is still robust in the proximity of the AFM re- 
gion, i.e., when there are large AFM fluctuations but no AFM 
order. What makes the mechanism robust is the fact that it 
is driven by confinement of topological defects (solitons) that 
are attached to each carrier (holes). In particular, as we will 
see below, this implies that the binding energy and the size of 
the pair are determined by different energy scales. 

We stait by deriving an extended Kondo lattice (KL) chain 
with a correlated (t — J) conduction band as the low-energy 
effective model H of the original two-band Hubbard Hamil- 



tonian H. The coiTelated nature of the conduction band is the 
main difference with the standard KL chain that was exten- 
sively studied in previous works |3|. With a simple analytical 
treatment of the fully anisotropic (Ising-like) limit of H, we 
show the origin of the two-hole bound state and the corre- 
sponding binding energy in the dilute limit. Our DMRG cal- 
culations allow us to extend these results to the fully isotropic 
(Heisenberg) limit and to finite hole concentrations. More- 
over, we show that the superconducting pair-pair correlations 
are the dominant ones over an extended and relevant region of 
the quantum phase diagram. Interestingly enough, the pairing 
remains robust in the absence of long-range AFM order (fully 
isotropic limit) indicating that a long enough AFM correlation 
length is sufficient for stabilizing the pairing mechanism that 
we discuss below. 

We will consider the following two-band Hubbard chain: 

where \ < j < L, L \s even, L + 1 = 1 [periodic boundary 
conditions (PBC)], rj = {l,u} denotes the lower and upper 
bands, and a — {t, i}. The diagonal energies are e„ = A„i/2 
and ei — —Aui/2 with A^,; > 0, and the density operators are 

From now on, we will assume that the mean number of 
electrons per unit cell is 1 < n < 2. For t^pf = and 
Uu > A„/, the ground state subspace, S, consists of states 
containing one electron per site in the lower band (only spin 
remains as a degree of freedom). In contrast, the sites of the 
upper band can be empty or singly occupied. In the strong 
coupling limit, A„; > tjj,y and Ui - A„; > tui, 

the low-energy spectrum of H is described by the effective 
Hamiltonian, H, that acts on the subspace S and results from 
applying degenerate perturbation theory to second order in the 
hopping terms: 

H — tuu 'y^.i'^l + lau^ia-u + ^^jau'^j+lau) + "y^.i^n ~ l>)njau 
jur) jv 
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where c 



j) (constraint of no double- 



occupancy), = {x,y,z}, = itl /U, 



^, and = 



"^tli/iUi - A„;) + 2tli/{Uu + Aui). Although the exchange 
interactions are isotropic ( = J,, and = Jk), we split 
the Heisenberg terms for reasons that will become clear later. 
Note that we have neglected the attractive, — ^fijUj^i, and 
the correlated hopping terms that also appear to second order 
in tuu to keep H simple and because they are not relevant for 
the pairing derived below. A more extensive study including 
the effect of these terms will be presented elsewhere 14]. To 
simplify the notation, we will introduce t — t^u, J — 
aJ = Jl ^ jy, Jh = Jf , and 13 J h = Jf = Jf ■ Where 
< a, /? < 1 determine the exchange anisotropics. H is an 
extension of the Kondo Lattice chain: the localized spins of 
the lower band are coupled via a Kondo interaction, Jk, to 
the itinerant electrons of the upper band which are described 
by a t — J model. In contrast to the usual Kondo Lattice, the 
itinerant electrons are strongly correlated and the local spins 
are AFM coupled via Jh- 

In the following, we set t = 1 as the energy scale and use 
J = 0.4, Jh ~ 0.5. These values correspond for instance to 
L/„ = 10, L/; = 16 and tu = \/2. We will start by assuming an 
Ising-like {(3 = 0) coupling between the localized spins to sta- 
bilize long-range antiferromagnetic order at T = 0. The fully 
isotropic case a = (3 — \ will be considered in the second 
part of the manuscript, so /3 = unless its value is explicitly 
specified. Below we will drop the band index in the definition 
of the different correlation functions because they are always 
applied to the conduction band. We will only compute the 
ground state properties. 

We use the DMRG method 15] |6l to study systems up to 
100 unit cells at T = 0. The PBC increase substantially the 
computational effort. In the finite-system step, we keep up to 
M = 1400 states per block and perform up to 12 sweeps. The 
weight of the discarded states is kept to the order 10^^ — 10^^° 
for /3 = 0, and 10"^ for /3 = 1 |7 1. 

We first consider the simplest case of Ising-like exchange 
interactions (a — 0) and a Kondo coupling much smaller than 
the rest of the terms in H. In this situation, the exchange Jh 
forces the localized spins to be AFM ordered in the ground 
state. For the Jk = ground state, each hole added to the 
conduction band carries a soliton or anti-phase domain wall 
(ADW) for the AFM order parameter ||8J. For Jk ^ 0, 
the single hole quasiparticle becomes a spinon-holon bound 
state [see Figjlja)] to avoid an energy increase of order LJk 
1^. The single added hole is topologically neutral: the holon 
and the spinon carry solitons with opposite "charge" (kink 
and anti-kink). We also note that unless the spinon and the 
holon are on the same site, the ferromagnetic link associated 
to the spinon increases the magnetic energy by J/2. This pro- 
vides an additional attractive force between the spinon and the 
holon. In this fully anisotropic limit, the spinon is immobile 
and the holon is localized around it ifTOl . 

The situation is qualitatively different for the two-hole 
ground state [see FigjTJb)]. The spinons attached to each 
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FIG. 1: (Color online) (a) Scheme of the single-hole quasiparti- 
cle. The holon and spinon carry solitons with opposite topological 
charges as indicated by the dashed line. The itinerant and localized 
spins are ferromagnetically aligned between the holon and spinon 
leading to a string or confining (linear) potential, (b) Scheme of 
the two-hole ground state. The holons carry solitons with opposite 
charge and are bound together by a magnetic string, (c) Formation of 
holon-holon bound state from two free holes. The two spinons with 
opposite spins and topological charges cancel each other, leaving two 
holons with opposite topological charges attached by the string, (d) 
Chsir) for a single hole. The arrow indicates the position of the 
spinon. (e) Chsir) for two holes with a — 1.0, Jk = 0.05, and 
L — 100. The reference holon is at r = 1. The arrow indicates 
the position of the second holon. Each holon carries an ADW as 
explained in the text. 

holon cancel each other (they have opposite spins and topolog- 
ical charges) leading to a bound state of two holons attached 
by a magnetic string. The cancellation of the two spinons low- 
ers the magnetic energy by ^ J. Consequently, we expect the 
two-holon bound state to be more stable than two independent 
spinon-holon pairs. This statement can be quantified by com- 
puting the exact ground state energies for one and two holes 
(this can be done for a = because the spins are not ex- 
changed by H). In particular, for Jk = 0.1, we get a binding 
energy As = Eg{Nh)+Eg{Nh-2)-2Eg{Nh-l) ~ -0.25 
for Nil — 2, where Eg{Nh) is the ground state energy for 
Nil holes. The holon-holon pair formation is illustrated in 
FigjTJc). The mutual cancellation of the two spinons leaves 
the two holons attached by a magnetic string. 

This simple picture for one and two holes has been dis- 
cussed previously in the context of a i — J model in a stag- 
gered magnetic field [11 1. In our case, the staggered mag- 
netic field h is not artificial because it is self-generated by the 
AFM ordering of the localized spins. As it was pointed out 
in ifTTI . the limit h (Jk ^ in our case) is singular: 
lim/i^^o^s = — J while the mean distance between the two 
holes diverges as h"^/'^ (Jk^^^)- This singular behavior is a 
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FIG. 2: (Color online) (a) Finite size scaling of As. The results in- 
dicate robust pairing in the thermodynamic limit for for all values of 
a. (b) Density-density correlation function Chh{r) for Jk ~ 0.1 
and a = 0.0,0.6,0.8, and 1.0. Chh clearly confirms the existence 
of a holon-holon bound state. Exact results are shown in open cir- 
cles for the Ising limit (a = 0) and the agreement with DMRG is 
excellent, (c) Chh(r) in the fully isotropic limit {a — (3 = 1) for 
Jk = 0.2, 0.3, and 0.35. A bound state is formed in this limit. 

manifestation of the qualitative difference between the single- 
and two-hole states: the binding energy Ab and the size of the 
two-holon bound state, Ip, are determined by two independent 
energy scales. While ^ — J for small enough h { Jk), Ip 
only depends on h (Jk) as long as J is nonzero. In particular, 
this shows that a negative binding energy, < 0, does not 
imply the formation of a two-hole bound state. The negative 
value of As for h ^ (Jk 0) just indicates that a single 
hole always creates a spinon [see Figjlja)] while this is not 
true for the two-hole state as shown in Figs[TJb) and (c). An 
infinitesimal field h (Jk) is enough for stabilizing the bound 
state due to the topological nature of the two-hole state: each 
hole carries a soliton and the two solitons become confined 
for any finite h. This remarkable property makes the pairing 
robust against the inclusion of a more realistic longer range 
Coulomb interaction in H. We will see below that this pairing 
mechanism survives in the absence of long range AFM order, 
i.e., when the effect of the localized spins cannot be replaced 
by a staggered mean field h because (S^,,) = 0. This is a very 
important qualitative difference relative to the case considered 
in 111 n. 

The above picture remains valid away from the Ising limit 
( a > 0). However, the hole can now move coherently in 
one sublattice. The process involves the exchange of the two 
spins on adjacent nearest and next-nearest neighbor sites of 
the hole via the a J term followed by two hoppings in the 
same direction. This coherent motion of the hole leaves the 
magnetic background unchanged. For < a < 1, the holon 
and spinon form a bound state. The magnetic structure of this 
quasiparticle is shown schematically in Fig. [ija). To confirm 
the above picture, we compute the hole-spin correlation func- 
tion Chsir) = {Sfn'l+iS^+r)^ where n'^ ^ 1 - {n^^ + n^i) 
is the hole density at site j. Figure[T|d) shows Chs{f) in the 
isotropic limit (a — \) with one hole in the conduction band. 
The correlator shows the spinon (indicated by the arrow) sep- 
arated from the holon (at r = 1) by a finite distance. Chs 
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FIG. 3: (Color online) (a) (a, u) phase diagram for L = 64 and 
Jk = 0.1. Near the Ising limit (a = 0.2), the pairing survives 
up to ~ 0.5 {Nh ~ 32). For a given q, the boundary between 
the two phases is determined from the condition As(a,!/c) = 
with Vc being the lowest value of v which satisfies this condition 
(holes are not bound for > i^c)- (b) As(z/) for a = 0.6. In the 
thermodynamic limit (L — + oo). As should be zero for v > u^, 
however As = 64) > due to finite size effects. 

shows clearly that the holon and spinon carry an ADW. When 
a second hole is added to the conduction band, the spinons 
cancel each other, as explained above, leaving the two holons 
attached to ADW's of opposite sign. This is also confirmed 
by the calculation of Chs [see Fig. [ije)]. The reference holon 
is at r = 1 while the average position of the second holon is 
indicated by the arrow (note that this average distance is artifi- 
cially increased due to the PBC). To the left side of the second 
holon, the spins at even numbered sites (r = 2 for example) 
are antiparallel to the reference spin; whereas to its right, they 
are parallel (e.g. r = 32). 

In the following we present numerical evidence of pairing 
and dominant superconducting correlations as a function of 
hole density v = Nh/ L. Figure|2|a) shows the pairing energy 



Ab = Eg{2) + Eg{Q) - 2Eg{l) in the dilute limit v 







versus 1/L and for different values of Jk and a. The finite 
size scaling of A^ shows robust pairing in the thermodynamic 
limit for all values of a. Figure[2]^b) shows the density-density 
correlation function Chhi^) = {^i'lT'i+r) foi" Jk = 0.1 and 
different values of a. Chh shows a clear maximum at a fi- 
nite distance confirming that the two holons are bound. As 
expected, the distance between the two holons increases with 
a. In the Ising limit, a = 0, the DMRG results are compared 
to the exact solution and the agreement is excellent. 

So far we have only considered the limit (3 — because 
the pairing mechanism is easier to identify in the presence of 
long-range AFM order It is natural to ask if the pairing sur- 
vives in the fully isotropic limit (a — (3 — 1) relevant for our 
original Hamiltonian H. In the absence of holes, the ground 
state of H only exhibits short-range AFM correlations due to 
the gap induced by the relevant Jk coupling. The Chh corre- 
lator for the two-hole ground state shown in Fig|2|c) provides 
clear evidence of the formation of a two-hole bound state. 
Moreover, the pair size Ip decreases monotonically with in- 
creasing Jk (slope of the confining potential) indicating that 
the pairing mechanism remains the same. 
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FIG. 4: (Color online) Single particle, Cm, and pair-pair, Csc, cor- 
relation functions. The circles show the numerical results, and the 
straight lines show the best fit. (a), (c) Jk = 0.05, a = 0.2. (b), (d) 
Jk = 0.05, Q = 0.6. CA/(r) shows an exponential decay whereas 
Csc{r) shows a power law decay. This confirms the emergence of 
dominant superconducting correlations. 

For finite hole concentration, a new length scale 1 /ly, that is 
the mean distance between holes, appears in the problem. We 
expect the previous analysis of the dilute limit to remain valid 
as long as /p ^ On the other hand, the pairing should 
be suppressed when the these two lengths become comparable 
because the effective interaction between holes is repulsive at 
short distances. This expectation is fully consistent with the 
numerical results shown in Figj3] The boundary between the 
paired and unpaired regions (see Figj3^) is shifted to lower 
hole concentrations when a gets closer to one, i.e., when Ip 
becomes bigger In particular, FigjSj) shows the evolution of 
the pairing energy As as a function of i/ for a = 0.6. Ab 
ceases to be negative for a hole concentration close to 20% 
{Vc ~ 0.2). This implies that if we vary the chemical potential, 
the states with odd number of holes are metastable as long 
as ly < Vc- We note that the critical hole concentration Vc 
remains significantly high {vc ^ 0.1) in the limit a = \. 

To confirm the presence of dominant superconducting 
correlations, we compare the single particle CM{r) — 

AJ) correla- 



'Eai^^r+rAa) ^d pair-pair Csc{r) 



tion functions, with AJ — -7^{cIj_t^c- 



y^\^i+iT^ll ^ Figure 
|4] shows results obtained for L = 80, Nh = 4, Jk = 0.05 
and a = 0.2,0.6. The single particle correlator decays ex- 
ponentially while the pair-pair correlator shows a slower al- 
gebraic decay. The exponential decay of the single particle 
coiTelator is due to the fact that holes are bound in pairs. In 
the dilute limit, the probability of finding the two holes sepa- 
rated by a distance r bigger than the pair size Ip (see Figpb) 
decreases exponentially in r. The algebraic decay of Csc(r) 
is expected for a Luttinger liquid of pairs. The comparison 
between Cm(?') and Csc{f) clearly indicates that the ground 
state has dominant superconducting correlations in the regime 
under consideration. 



A few comments are in order. We verified that the pair- 
ing (Ab < 0) persists for smaller values of Jh such as 
Jh — Jk = 0.1 . The small Jh regime of H is relevant for 
describing lanthanide and actinide based compounds in which 
localized f-electrons interact via Kondo exchange with elec- 
trons in the conduction band. Conduction bands with strong 
Sd-character are correlated and would provide a natural real- 
ization of our H. If the t—J conduction band of H is replaced 
by the original Hubbard upper band of H, one can study the 
evolution of our pairing mechanism as a function of [/„/t (the 
standard KL model is recovered for Uu/t = 0) 14 1. 

In contrast to pairing mechanisms driven by an attractive 
short range interaction (like the AFM exchange J), our mech- 
anism is robust under the inclusion of a more realistic longer 
range Coulomb interaction. This results from the fact that A^ 
and Ip are determined by two independent energy scales (J 
and Jk)- In other words, a big enough value Ip reduces the 
effect of longer range Coulomb terms without reducing the 
value of As (note that Ip and As are anticoiTelated when the 
pairing is produced by a short range attractive potential). 

Finally, our pairing mechanism should also persist for 
weakly coupled chains (small inter-chain hopping) due to its 
topological nature. While single-particle coherent inter-chain 
hopping is not possible due to the soliton that is attached to 
each hole, a pair can hop coherently between chains because 
it is topologically neutral (soliton-antisoliton bound state). 
Such coherent pair-hopping should stabilize a superconduct- 
ing state below a finite critical temperature. A detailed study 
of this extension to higher-dimensional systems will be pre- 
sented elsewhere [41. 
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